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O ' I- INTRODUCTION 

o ■ 

, The derivation of the Casimir effect in its original context iJj] was the calculation of the electromagnetic force 
between parallel perfectly conducting plates placed in a vacuum. The geometry of a cavity formed by the parallel 
^ \ plates is also the best context to manifest the nonstationary (dynamical) Casimir effect. The theoretical calculations 
'^<^ • were first carried out by Moore The force between plates is dynamically modified, but one encounters more 
\^ \ amazing phenomenon being the vacuum emission of photons with nonclassical properties. The latter topic became 
the subject of intensive studies since it was recognized that the production of photons could be significantly enhanced 
under the parametric resonance conditions Fr om the experimental point of view, the hopes for observing 

vacuum radiation are based on the belief that vibrating cavities can be adjusted to parametric resonance conditions. 

Quantum parametric resonance is a complex process and traditional and well-tried methods for calculations in 
quantum electrodynamics are obsolete here. We must examine open quantum systems under the action of external 
periodic perturbations. Moreover, the focus of our interest are the instable systems. If we describe physical quantities 
I in resonant systems, then the corresponding mathematical functions inevitably yield narrow peaks or steep stairs. 
. A naive approach to make perturbations in small amplitudes of cavity vibrations breaks down. Therefore, we must 
' go beyond the standard mathematical treatment known from QFT and pursue some new machinery. Most ways 
of tackling the problem follow the quantization scheme described by Moore |2] in case of one-dimensional cavities. 
^ ^ Then, special mathematical techniques are used to trace the evolution of the quantum system with time-dependent 
boundary conditions. Dodonov and his collaborators reduced the analysis of field equations in a cavity system to 
■ some infinite set of coupled and time-dependent oscillators. Another idea based on a concept of renormalization group 
\ was proposed in 0] . Our suggestion Q is to trace the evolution of the quantum fields in a vibrating cavity system 
^ . following classical optical paths. It is related to the iteration procedure described in and techniques developed 
^ '' for classical field theory in and In fact, the resonant evolution of the physical system is similar in classical 
^ . and quantum models. Recently, a similar idea to g round calculations on the analysis of classical optical paths was 

advocated for a stationary Casimir effect in p^^J I ^ . 
^ . , In the present paper, we apply our technique |g| to examine cavities with various kinds of their motion and different 
' types of boundary conditions. The calculations are involved in determining the essential physical outcome. In 
particular, the conditions necessary to keep the parametric resonance and the critical exponents describing resonant 
evolution are derived. All results are in agreement with those obtained for cavities with mowing walls using other 
techniques. We give more information about resonance windows, more subtle structure of energy densities when 
amplitudes increase or the adjustment of resonance parameters is getting worse and other new details. 
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II. A CAVITY WITH ONE OSCILLATING WALL 



Most researches on vibrating cavities concentrated on a model of a cavity formed by a static wall and a moving one. 
The static wall is fixed at x = 0, while the moving one follows its trajectory x — L{t). The condition that L{t) > is 
assumed. It is convenient to make use of the billiard function 14]: 

f{t + L{t))=t-L{t) . (1) 
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This function contains all information about possible scatterings of massless particles from the moving wall. It can 
be also understood as the trajectory of the wall in light-cone coordinates. Its derivatives are recognized as retarded 
Doppler factors (the retardation relation is t = t* + L{t*)): 

m ^ '-Mn . (2) 

Consider now a classical scalar field A(t, x) in a one-dimensional cavity. It obeys the wave equation: 

{d^,^dl)A{t,x)^Q. (3) 

The field is subject to some boundary conditions. Usually, we assume either Dirichlet's boundary conditions: 

A{t,G)^A(t,L(t))=0, (4) 

or Neumann's boundary conditions at the instantaneous positions of walls: 

d^Ait, 0) - {L{t)dt ~ d,)A{t, L{t)) = 0, (5) 

In both cases, the classical solutions are given by: 

A{t,x)=ip{t + x)T^{t~x), (6) 

where the upper sign corresponds to Dirichlet's case, while the lower one to Neumann's case. The profile function is 
subject to the following relation: 

V{r) = (^(/(r)) . (7) 
The energy density of the classical wave packet and its total energy are given by: 

Too(i, x) = l/2{dtAf + l/2{d.,Af = g{t + x) + g{t - x) , (8) 



.L(t) nt+Lit) 

E[t)^ dxToo{t,x)^ dr g{T) , (9) 

JQ Jt-L(t) 

where q{t) = </5^(t). Let us now describe classical optical paths. We consider a massless particle that starts at time 
T moving right from the position x = Q. Denote by T„(t) and T*(t) times of n-th collisions of the particle with the 
moving wall and the static one respectively. It is easy to notice that r„(T) = where the notation (/~^)°"(t) 

means here and throughout this paper n-fold composition o o . . . o /~^, while T*{t) are to be determined from 
the retardation relations T*{t) -\- L{T*{t)) = Tnir). The parametric resonance is related to the existence of periodic 
particle trajectories 8, llj. Each periodic trajectory obeys a following condition (for any non-negative integer n): 

Tn{To) = TO + nT , (10) 

where tq is a starting point in time and T is a period. The retarded times are respectively T*{tq) — T„(to) — T/2. 
Any periodic particle trajectory appears on condition that the mirror trajectory has return points: 

L{t; + nT) = T/2 . (11) 

It is usually assumed that T/2 refers to the length of the static cavity, so that T/2 = L = L{0). The crucial 
characteristic of a particle trajectory is the cumulative Doppler factor Q: 

The above function has a direct physical interpretation. If a massless particle starts at time r moving right inside 
the cavity, then after n rebounds from both cavity walls its kinetic energy changes by the factor Z?„(r). A trajectory 
is called positive (or negative respectively) if its cumulative Doppler factor tends to infinity (or zero) with increasing 
n. We denote by t+ and r_ starting points for positive and negative periodic trajectories. If positive trajectories 
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exist, then any trajectory approaches one of the positive trajectories for long times. The same is true for negative 

trajectories if we consider the evolution backward in time. In other words, T'„(t) Tn{T+) and T~^{t) — > T~^(t_) 
for large n. 

The evolution of the energy density can be deduced from the following formula, 

f?(T„(r)) = g{T)Dl{T) , (13) 
together with the asymptotic approximation for long times (i.e. n ^ 1)0: 

ginT + T++ eD-\T+)) - ^?(r+ + e)Dl{T+ + e) , (14) 



QinT + T_ + e) ^ f?(r_ + ei^„(T_))i?2 (r_ + eI?„(T_)) . (15) 



The above formulas explain the formation of travelling narrow packets in the energy density Too(t,a:) |3l I5lll5| . The 
profile function q{t) has peaks located at the spots of positive periodic trajectories t+. The height of a peak grows 
like r'^(T+), while its width diminishes like Z),7^(t+). The value of the energy density near the peak decreases like 
Z?^(t_). The evolution of the total energy can be analyzed using the formula (we assume here tq < L): 

rTn(To) pro 
E{n{To)) = / dT b{t) = / dT B{T)Dr,{T) . (16) 

JT^-iira) Jf(ro) 

To calculate the total accumulated energy inside the cavity after some time we need only an initial energy density 
(i.e. gi(r) for t < L) and the knowledge of cumulative Doppler factors. In fact, it is usually enough to know these 
factors in the vicinities of spots of periodic trajectories. 

Now, we turn into the basic model of vibrating cavities specified by the sinusoidal cavity wall motion: 

L{t) = L + ALsm{ujt) , (17) 

where AL < L and luAL < 1. The parametric resonance frequencies are lun — Nn/L, where N is the order of the 
resonance. Positive and negative periodic particle trajectories correspond to the following starting points: 

{-N + 2m + l)L {-N + 2m)L . , 

•^+"1 = , T_,„ = — , (18) 

where m = 0, 1, A'' — 1. The cumulative Doppler factors can be calculated: 

Dn[T±m) = -; -ry 

= exp {±2nujNAL) for AL < i . (19) 

We conclude that there are N travelling peaks in the energy density. The average distances between peaks are equal. 
For N even, N/2 peaks are moving right and N/2 peaks are moving left. For odd N the number of peaks of each 
kind differs by one. All the peaks have the same height and the same width. The total energy grows exponentially 
with time. These classical results agree with the results of the corresponding quantized cavity model 0. This fact 
will be explained later. 

The above picture is true provided that the amplitude of cavity oscillations is small. Strictly speaking, we need 
AL/L < N. For larger amplitudes, the behavior of the cavity system is slightly modified. If AL/L > M/N 
{M = 1, 2, 3, ...), then there appear additional positive trajectories corresponding to the following return points ((TT|l : 

WAfALsin (watt^^) = Mtt , (20) 

and the corresponding cumulative Doppler factors are: 



It follows that beside the first series of N peaks in the energy density, there appear next M series of N peaks each. 
The peaks of additional series are smaller and located between higher peaks. 
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For practical purposes, it is important to consider off resonant behavior of the cavity system. The instabihty in 
the field theory appears not only for finely tuned frequencies. We are usually dealing with some band structure. To 
establish the resonance window, we move the resonant frequency: oj = ojn + Aw. The equation for return points Hll|l 
is now: 

. , LAuj 

sin wr* = — , 22 

ujAL 

so that solutions exist provided that: 

Aw AL 

— < — . (23) 
uj L 

This condition agrees with the result of the quantized model [l6j | . It defines the windows for the resonant frequencies 
to keep the resonant instability. The cumulative Doppler factors at starting points of periodic trajectories are: 



Again, the above picture is exact provided that AL/L + (1 + 1/N)Alu/u! > 1/N. Otherwise, next series of smaller 
peaks in the energy density will appear. 

Now, we present quantum formulas. The vacuum expectation value of the energy density is given by 0: 

{Too{t,x)) ^ git + x) + git~x) , (25) 

where 

^'M = -i^^'(^)-i'^t^Kr), (26) 
and the Schwartz derivative is defined by: 

S[«l^|-|(f)". (27, 

The function R{t) is subject to Moore's equation: 

RiT) - i?(/(r)) = 2 , (28) 

where our billiard function appears in the quantum case as well. The formula for the evolution of the energy density 
is given by: 

^?(T„(r)) = g{T)Dl{T) + A,,{t)DI{t) , (29) 

where the cumulative Doppler factor is defined again by the same formula Eq. (|12|l and the definition of the quantum 
term called the cumulative conformal anomaly contribution Q is given below: 

1 1 " 

^"^^^ = 24^^[^"](^^ ^ -^E^fc"'W'5[/](rfc(r)) . (30) 

fc=i 

The quantum formula for the total energy accumulated in the cavity is: 

/■T,.(ro) 

E{n{To)) = / dT g{T) 

dr [g{T)+A,,{T)]Dn{T) . (31) 

/(to) 

The quantum anomaly implies the difference between classical and quantum formulas for the evolution of the energy. 
However, under the resonance conditions the leading role is played by the cumulative Doppler factor. This fact 
explains why most of classical results match quantum ones. It refers to the band structure of resonance frequencies, 
the formation and the shape of travelling packets in the energy density, the exponential growth of the total energy. 
To make predictions for these features, the classical theory could be good enough. But it is not true for the lowest 
resonance channel. In this case, the anomalous mechanism of energy growth clashes with the resonant enhancement 
of the initial vacuum fluctuations . The lack of the instability for the lowest resonance frequency is also related to 
the SL{2, R) symmetry of the cavity system which is the remainder of the conformal group. 

Finally, let us also conclude from the above results that there is no difference between imposing Dirichlet's and 
Neumann's boundary conditions in the context of our cavity models. 
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III. A CAVITY WITH TWO OSCILLATING WALLS 

A natural generalization of the considerations of the previous section is to allow both walls to oscillate 0, 0, ^3 ■ 
Consider now that the right cavity wall moves along the trajectory x ~ Li{t), while the left one follows x = —L2{t). 
We assume also that the cavity was static in the past, so that Li{t) = L2{t) = L/2 for t < 0, and L corresponds to the 
static cavity length. The general formulas given in this section will be rigorously justified if we assure that Li{t) > 
and L2{t) > (the walls never collides) and the wall velocities are never close to the velocity of light. The billiard 
functions are implemented for both trajectories (i = 1,2): 

Mt + L,{t)) ^ t - L,{t) . (32) 

Moreover, it is convenient to define an additional pair of billiard functions: 

h = /2 ° /l , fR = /l O /2 . (33) 

They describe double reflections of massless particles from the cavity walls. The set of classical optical paths is now 
more complex. Consider again a massless particle that starts at time t moving from the position x = 0. Denote 
by TLn(T) (and Tjinir) respectively) times when the particle returns to its initial position a; = after n rebounds 
from both walls assuming that the first collision is with the left (or with the right) wall. It is easy to find that 
T£(fl)„(r) = (/L(ii) )""(''')• Further, T'2(i?)n(''') ^'^'^ -^L(fl)n(''") times of successive collisions with the left (right) and 
the right (left) walls for both types of trajectories respectively. They can be derived from the following relations: 

TUR)u{r) + ii(2)(T2(«)„(r)) - TLiK)jT) , (34) 



nUnir) + i2(l)(n(fl)n(^)) = ./i(2)(Tl(«)„(t)) . (35) 

For both types of trajectories cumulative Doppler factors are defined: 

1 " 

DL(R)n{T) = -^r, — = TT iL(R){TL(R){T)) 

J-L(R)n{T) 



^ l-ii(2)(7^LV).W) A l-^2(i)(T2f;,),(r)) 



T-r -n.^n^L(fl)fcv y-r ""^U 

M i+w^w^ L 1 1+^2(1 



(36) 



Classical solutions of a wave equation inside a cavity with two moving walls are given by: 

A{t,x)^ipL{t + x)TVR{t-x), (37) 

where again the upper sign corresponds to Dirichlct's boundary conditions and the lower one to Neumann's boundary 
conditions on the walls. The profile functions in both cases obey the same relations: 

Vl{t) = ^fl(/i(T)) , ^fl(T) = ^L(/2(r)) . (38) 

The energy density of the classical wave packet and its total energy can be represented as: 

Too{t,x) ^ QL{t + x) + QR{t-x) , QL(K){'r) ^ ^ (39) 

j-Li(t) pt+Li{t) rt+L2{t) 

E{t)= dxTon{t,x)^ dTgL{T)+ dT qr{t) . (40) 

J-L2{t) Jt-L2(t) Jt-Li(t) 

The classical evolution of the system can be traced with the help of the following formulas: 

QL(R){TL(K)n{T)) = Qur){t)DI^„;^^{t) , (41) 



E{t) = f dT QL{T)DLn{T) + f dr qr{t)Dr,,{t) , 

Jh{r2) Jfl(ri) 



(42) 
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where ri(2) = + ^i(2)(i))- If we take n big enough to support that ri(2) < L/2, then we obtain the formula 

that allows us to calculate the evolution of the total energy from the knowledge of Doppler factors and the energy 
distribution in the initial state. 

Before we turn into examples of cavity systems, we explain why many specific systems of cavities with two oscillating 
walls just reproduce results similar to those obtained in cavity models with one oscillating wall. The reason is that 
billiard functions (|33|l correspond to some physical trajectories in two dimensions. It means that we can find a 
trajectory x — L(t) which describes some realistic motion of the wall such that: 

fLiB)(t + m) = t-L(t) . (43) 

We can verify that L(i) > 0, L(i) = i for t < and L(t) is never close to the speed of light. Obviously, in general 
such "effective" trajectories are different for and /j^, but in both cases the following composition rules are satisfied: 

L(t) = + U{t2) , (44) 

l-L(t) l-ii(ti) 1-1.2(^2) 



1 + L(<) 1 + Li(ii) l + i2(t2) 



(45) 



where t ± Lit) = ti ± L\{t\) and t\ =F L\{t\) = <2 ± -^2(^2) for /l (the upper sign) and f^. respectively. 

The latter relation Ea. l|45(l states just that L(t) is the relativistic sum of velocities L\{t\) and L^if-'i)- However, this 
composition rule has nothing to do with any relative motion since there exists no inertial frame where we could make 
respective time moments synchronous. Nevertheless, both left- and right-movers in cavities with two oscillating walls 
correspond to some models of cavities with a single oscillating wall. In other words, each mode of a cavity system 
with two moving walls can be described by a simpler model discussed in the previous section. 

Let us discuss specific types of symmetric and antisymmetric oscillations of cavities. In the paper of Dalvit and 
Mazzitelli 0, these types of cavity oscillations are called breathing modes ("electromagnetic antishaker") and trans- 
lational modes ("electromagnetic shaker"). For symmetric oscillations we have L\(t) — L^it) and the cavity oscillates 
symmetrically with respect to the center of the cavity. The billiard functions and Doppler factors are identical, namely 
/i = /2, jh = Ir and Dl = D^. The model is equivalent to the cavity model with one oscillating wall, where the 
prescribed trajectory of the wall is defined by the billiard function /i o fi. For antisymmetric oscillations we have 
Liit) + L2{t) = L, and the cavity is oscillating as a whole with its length kept constant in the laboratory frame. 
Then, /i(2)(t) = f2{i)i''' — L) — L and /l(_r,)(t) = /^2)(/i(2)(''') ~ L) — L. It is interesting to note that for periodic 
translational cavity oscillations with L being a period, the billiard functions Jl and fa correspond to static ones. It 
follows that there is no effect of the cavity movement on the energy distribution of fields and the whole cavity system 
is equivalent to the static one. 

Finally, we analyze harmonic oscillations of cavity walls: 

Li{t) = I + ALi sin (tjfli) , 
L2{t) = ^ + AL2 sin {ujLt -5) + AL2 sin (5 . (46) 

Again, we want to establish resonance windows, find the evolution of the total energy and describe the shape of the 
local energy density. We must first detect positive and negative periodic trajectories of massless particles moving 
inside the cavity. They can be derived from the generalization of Eq. Hll|) : 

Li(ii) + L2(ti±^)-| , (47) 

where the upper (lower) sign corresponds to the "left" T^n ("right" Tr„) trajectories. The cumulative Doppler factors 
at the starting points of periodic trajectories can be calculated then: 

1 — ujjjALi cos {ujRti) \ " / 1 ^ ujlAL2 cos (wl^i T ojlT/2 — S) 



DL(R)n yi +^^ALiCOs(w_Rtl)y \\+UJLAL2COs{uJLtlTi^LTl2- S) ^ 

= exp {-2n[ujn,ALi cos (w^ti) + ujlAL2 cos (w^ii T i^LT/2 - S)]) for AL; < L . (48) 

First, consider resonant oscillations with lol = i^r = = Nt: / L. If there is no dephasing (5 = 0, the system is 
equivalent to the cavity with one oscillating wall. The oscillation frequency is w^r, the effective cavity length is L and 
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the effective amplitude of oscillations is ALi + (— 1)^AL2- For ALi = AL2 and TV odd the effective cavity system is 
static. The same is true for ALi = AL2, S = n and N even. These results agree with the paper 0. If the amplitudes 
of oscillations are equal ALi = AL2 = AL and S 0,tt, then the solutions of Ea. (|47|l are ti — {SL/tt + 2k)L/N and 
ii = {N — 2k ~ 1)L/N, where fc = 0, 1, iV — 1. The respective cumulative Doppler factors are: 

D.n{Tk±) G^^exp(±2nwArAL(l + (-l)^cosJ)) . (49) 

There are N positive and N negative particle trajectories. It follows that there are N peaks in the energy density. 
The calculations for ALi ^ AL2 are a bit more complicated. However, if ALi + (— 1)^AL2 cos (5 7^ then there exist 
N positive and N negative particle trajectories. The corresponding cumulative Doppler factors are given by: 

Dn = exp {±2nujN{ALi + (-1)^AL2 cos^)) . (50) 

For ALi + (— 1)^AL2 cos 5 = there is no exponential instability in the cavity system. 

To establish the resonance window we assume = u)r = cu^ + Aw. The relation Ea. H47|l yields: 

ALi sin (wii) + (-l)^AL2sin(wti - (5) + AL2 sin (5 = . (51) 

id 

To keep the parametric resonance conditions for the vibrating cavity Eq.l^HJ it is enough to support that: 



Aiu AL2 smS+ V(A£i)2 + (A£2)" + 2(-l)^A£iA£2 cos<5 

~ < L ■ ^^^> 

The quantum formulas for cavities with two oscillating walls include anomaly contributions in the respective way. 
Since the modifications are straightforward, we skip their discussion here. 



IV. CONCLUSIONS 



An optical approach to the dynamical Casimir effect based on the construction of classical optical paths is useful. 
In particular, it is a powerful technique to analyze resonance systems. We can trace the behavior of a cavity system 
with moving walls when we increase amplitudes or loose fine tuning for resonance eigenfrequencies. We can calculate 
modified critical exponents, find more subtle structure of travelling peaks in energy densities and control the conditions 
for the instability of the system. We have established some "composition laws" coming into effect if we want a more 
complicated model of a cavity with two moving walls reduce to the analysis of a cavity with one moving wall. Moreover, 
our understanding of the mechanism of parametric resonance in cavities with moving walls and the implications of 
quantum field theory is much better. For example, we can judge when some semi-classical approximation is accurate 
or we can explain why some resonances are not observed at the quantum level. It is further remarkable that we can 
distinguish between exponential and power-like instability and observe the role of symmetries at the quantum level 

We have seen in several examples that many qualitative or even quantitative conclusions can be drawn from classical 
versions of respective cavity models. Presumably, the parametric resonance mechanism and evolution laws are similar 
in both cases. However, we should be careful about quantum anomalies that can significantly modify our predictions. 
A subtle role of quantum modifications can be traced within our approach. 

The results derived for Dirichlet's and Neumann's boundary conditions for one-dimensional cavities are the same. 
This point was already made in The cavities with two oscillating walls possess a double set of characteristic 

functions, namely two families of both cumulative Doppler factors and cumulative anomaly contributions. It is related 
to the fact that from a given space-time point there are two independent directions to follow trajectories of massless 
particles. In three dimensions, such functions depend also on localizations and orientations. 
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